Centrality dependence of charged hadron and strange hadron elliptic flow 
from ,/s7r = 200 GeV Au + Au collisions 
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We present STAR results on the elliptic flow V2 of charged hadrons, strange and multi-strange 
particles from y^s^ = 200 GeV Au + Au collisions at RHIC. The detailed study of the centrality 
dependence of 172 over a broad transverse momentum range is presented. Comparison of different 
analysis methods are made in order to estimate systematic uncertainties. In order to discuss the 
non-flow effect, we have performed the first analysis of V2 with the Lee- Yang Zero method for Kg 
and A . 

In the relatively low pr region, pr < 2 GeV /c, a scaling with mr — m is observed for identified 
hadrons in each centrality bin studied. However, we do not observe V2{pt) scaled by the participant 
eccentricity to be independent of centrality. At higher pr , 2 GeV/c < pr < 6 GeV/c, 112 scales 
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with quark number for all hadrons studied. For the multi-strange hadron Q, which does not suffer 
appreciable hadronic interactions, the values of V2 axe consistent with both mr — m scaling at low 
Pt and number-of-quark scaling at intermediate pt- As a function of collision centrality, an increase 
of pT-integrated V2 scaled by the participant eccentricity has been observed, indicating a stronger 
collective flow in more central Au + Au collisions. 

PACS numbers: 25.75.Ld, 25.75.Dw 

I. INTRODUCTION 

The event azimuthal anisotropy with respect to the reaction plane has been widely studied in order to evaluate 
the collective behavior of the matter produced in high-energy nuclear collisions P, 0, [1] . The initial configuration 
space anisotropy is expected to be sclf-qucnched by expansion and reduced by frequent re-scatterings in the hot 
and dense medium created in such collisions. The final observed momentum space anisotropics, therefore, carry 
information about the early stage collision dynamics [1, @|. The experimental results of the second harmonic 
azimuthal anisotrop y, e lliptic flow, V2^ from Au + Au collisions have demonstrated the development of partonic 
collectivity 0, H, [a. Il0l|. Further detailed analyses of the hadron mass dependence of V2 suggest that the system 
has been in the deconfined state with constituent quark degrees of freedom prior to hadronization [ll|, [l3| . 
Furthermore, results of multi-strange hadron transverse momentum distributions and V2 indicate that the system 
reached thermalization at the partonic stage [3, [l^, [l^, [13] ■ 

Hydrodynamic model calculations, with the assumption of ideal fluid behavior (no viscosity), have been successful 
when compared with the experimental data at RHIC [1, H, [HI, [H] ■ It should be noted that up to now, the discussions 
of the underlying dynamics of the thermalization at RHIC are inconclusive. Some initial evidence for thermalization 
was provided by the quantitative agreement of V2 results between ideal hydrodynamic model calculations and data for 
identified hadrons tt, K, p and A [ll|, [l^ from minimum bias (0 - 80% centrality) Au + Au collisions shown 
in Rcfs. [3, [^, ideal hydrodynamic model calculations have failed to reproduce the centrality dependence of tt and p 
V2 in Au -f Au collisions. In addition, the discussion based on the integrated y2/£vart of charged hadrons suggests 
possible thermalization only for the most central collisions at RHIC (see Refs. [2l|, [23] and references therein). Here 
the participant eccentricity, e parti is the initial configuration space eccentricity of the participants. From peripheral 
to the most central Au + Au collisions, the values of V2l£part increase as a function of the scaled charged hadron 
multiplicity, as predicted by a model calculation assuming the low density limit [23j of single forward nucleon-nucleon 
collisions. This analysis indicates that the system has probably not reached thermalization for most peripheral Au+Au 
collisions. 

Hydrodynamic model calculations predict a characteristic dependence of the elliptic fiow and transverse momentum 
spectra on particle mass and collision centrality. Nevertheless the comparisons made so far have been mostly restricted 
to identified hadrons from minimum bias collisions or integrated V2 of charged hadrons [2l| . Systematic comparisons 
for identified hadrons at different collision centralities are still scarce ^ . In order to fill this gap and further advance 
our understanding of the properties of the medium created in high-energy nuclear collisions, in this article we report 
the centrality dependence of the azimuthal anisotropy parameter V2 (elliptic flow) in ^/s^ ~ 200 GeV Au + Au 
collisions. The centrality dependence of V2 for identified hadrons A'^, A, S, and fl, and the scaling properties as a 
function of number of quarks within a given hadron and the transverse kinetic energy my — m are reported. Results 
from the Lee- Yang Zero method [13, HH] for unidentified charged hadrons {h^ + h~), and A are also reported. 
In complex nuclear collisions different systematic errors on V2 can arise from different analysis methods. In this 
paper, the systematic errors are analyzed b y co mparing the standard event plane method [3, [M] with results from the 
Lee- Yang Zero [l^, four-particle cumulant [2^. [26t. and 77 subevent ^ methods. 

The paper is organized in the following way: we discuss experimental cuts, data selections and methods used for 
unidentified charged hadrons and identified hadrons in Section [TT] Section IIIII gives the results on V2 for unidentified 
charged hadrons and identified hadrons along with a discussion of the systematic errors extracted from different 
analysis methods. Section HVl presents a comparison with model calculations as well as a discussion of scaling and the 
systematics of the V2{pt) distributions from SPS (y^s^= 17.3 GcV) and RHIC (^^8^= 62.4 and 200 GeV). Finally, 
the summary of the analysis and the outlook are presented in Section [Vj 

II. METHODS AND ANALYSIS 

In this paper, if V2 is used without curly brackets it is an abbreviation for W2{£'-P2}, that is, V2 relative to the second 
harmonic event plane. The systematic uncertainty from the event plane resolution is constant at each centrality and 
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expected to be smaller than that of the observed differential flow, and is not folded in. Other systematic errors are 
discussed in Section Hill 



A. Data sets 



For this study, the high statistics data from ^^s^ = 200 GeV Au + Au collisions collected by the STAR experiment 
during RHIC's fourth year (2004) of data taking were analyzed. STAR's main time projection chamber (TPC) [l^l 
was used for tracking and identification of charged particles. The TPC records the hits used for reconstructions of 
the tracks of the particles, enabling the measurement of the momenta of the particles, and identifying the particles 
by measuring their ionization energy loss. The TPC provides complete azimuthal coverage and complete tracking for 
charged particles within ±1.8 units of pseudorapidity. The two forward time projection chambers (FTPC) cover two 
sides of the collision with 2.5 < | | < 4.0. The FTPCs also provide a tool for studying non-flow effects. About 
25 million minimum bias events (0-80% most central of the hadronic interaction cross section) were analyzed, which 
increased the statistics for flow analysis by more than a factor of 10 compared to the previous measurements [lol. [T3| . 

The centrality definition of an event was based on the number of charged tracks in the TPC with track quality cuts, 
which are |?7| < 0.5, a distance of closest approach to the primary vertex (DC A) less than 3 cm, and fit points more 
than 15. These events were grouped into three centrality bins, which were central (0-10%), mid-central (10-40%), 
and peripheral (40-80%). In addition, the central dataset (0-10%) was enhanced by online triggering on the most 
central events with the Zero-Degree Calorimeter (ZDC) [l^, thereby getting an additional ^^19 million events for a 
similar centrality bin. Within the statistical uncertainty, the results from the central trigger dataset were consistent 
with those from the minimum bias trigger. 

The analyzed charged particles were identified as the track helix in the TPC magnetic field. The charged tracks 
were selected with a transverse momentum range of 0.15 < pt < 2.0 GeV /c unless indicated otherwise, and a pseu- 
dorapidity range of | ?7 | < 1.0. For the Lee- Yang Zero product generating function analysis the r] interval increased to 
\ 7] \ < 1.3 in order to obtain more particles. A minimum of 15 track fit hits and a ratio of hits to maximum possible 
hits > 0.52 was also required. To improve the selection of good tracks from the primary collisions, the distance of 
closest approach of the analyzed tracks to the event vertex had to be less than 2 cm. Tracks of charged daughter 
particles stemming from weak decay, which tend to be at large distances, are not subject to this cut. 



B. Particle identification 

We identified , A (A), S~ (S"*"), (11''') through their decay channel: 7r+ + 7r~, A p+n" (A p+n+), 

E- ^ A + n- (S ^ A+ TT+), Q- ^ A + K- K+). The charged pions, kaons and protons were 

identified via their energy loss in the TPC According to the (multi)strange particle decay properties, topological 
cuts and kinematic cuts were applied to reduce the combinatorial background. The detailed description of the analysis 
method can be found in Refs jT^, [2§| for and A , and Refs. [H, [1^ for S and SI. 

Figure [T] shows the invariant mass distributions for (al) , (bl) A -I- A , (cl) + and (dl) Sl^ + fl^ 
for a given px bin from ^/s^ = 200 GeV minimum bias (0-80%) Au + Au collisions. Clear signal peaks are seen 
at the values expected for the particle mass above the combinatorial background. The measured invariant mass 
distributions contain both signal (Sig) and combinatorial background {Bg). For Kg and A , the measured invariant 
mass distributions were fitted by a polynomial (up to 4*'* order), which represents the background, and a double- 
Gaussian function, which represents the signal. The double-Gaussian was used because of tails on the distribution. 
For multi-strange baryons S and SI, the Bg was estimated by rotating the transverse momentum of the daughter A by 
180°. This operation breaks the correlation between the A and the other daughter particle. The resulting invariant 
mass distributions provide a goo d ap proximation for the true background distribution. The detailed description of 
the method can be found in Ref. [12 |. 

For V2 of the identified particles, , A -I- A , -|- S and Sl^ + Vl , the V2 versus niinv method is used in this 
analysis 31|. Since V2 is additive, one can write the total V2^^^^^ as a sum of Sig and Bg contributions weighted by 



their relative yields: 



wf'^''''^^(m„i,) = uf*^ • — (TOi„„) + V2^{miny) ■ -—^^——{m,ny). (1) 

big + Bg big + Bg 



This method involves the calculation of V2^^'^^^ as a function of TOm„ and then fitting the distribution using Eq. ([T|) 



with measured relative yields and parameterizations of u^*" and v"^ {rriinv)- The „ . (niinv ) distribution is the 
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FIG. 1: (Color online) Plots (al)-(dl) represent the invariant mass distributions for Kg (1.4 < Pt < 1-6 GeV/c), A (1.4 < pt < 
1.6 GeV/c), H (2.3 < pt < 2.6 GeV/c), and (2.5 < pr < 3.0 GeV/c), respectively, from ys^= 200 GeV minimum bias 
(0-80%) Au-\-Au collisions. The dashed lines are the background distributions. The corresponding data for the V2 distributions 
are shown in plots (a2)-(d2) as open circles. The thick-dashed, thin-dashed and the dot-dashed lines represent the relative 
contributions of V2(Sig), V2{Bg), and V2{Sig + Bg), respectively. For clarity, the invariant mass plots for Kg, A, H, and Q,, are 
scaled by 1/50000, 1/170000, 1/2.5, and 1/3, respectively. The error bars are shown only for the statistical uncertainties. 



Bg divided by {Sig + Bg). The g^^'^^g {'minv ) distribution is simply calculated by 1 — si^+Bg The term 

{rriiny) is parameterized as a linear function in order to take care of the non-constant value as a function of 

rriinv The fit result ff*^ is the final observed V2- The fit results for Kg , A -I- A , 5~ -I- 'E^ and VL~ -\- VL^ are shovifn, 
as dot-dashed lines, in Fig. [T](a2), (b2), (c2), and (d2), respectively. Note that the anisotropy varies as a function 
of Pt and hadron mass. In this figure, the V2 are shovirn for different hadrons with the px cuts that are listed in the 
caption. How this method works well for measuring signal V2 is explained as following: a set of data points are used 
in the fit over a wide mi„„ region for Sig and Bg. Data points far from the mass peak constrain (niinv) since pure 
Bg is expected in this region. (The disagreement at M,;„„ ~ 1.29GeV/c^ in (cl) is caused by mis-identified hadrons 
in the Lambda-pion combinations, which is explained in Ref. [l^.) Under the peak, the wf'^^^^(TOmi') is dominated 
by the Sig distribution. Finally, the V2 signal is extracted by the fitting method shown in Eq. [T] 

The results obtained with this technique are in good agreement with the ones from the method used previously 0] . 
Note that the subtraction procedure used to extract the V2 signal for a given identified particle is independent of the 
flow correlations. The V2 distributions of the overall signal and background are evaluated by a specific flow analysis 
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method. These methods will be discussed in Section fll CI 

C. Flow analysis methods 

The systematic uncertainty of the Event Plane method is evaluated by comparing the results to those obtained by 
other techniques for measuring anisotropic flow. The various methods have different sensitivities to non-flow effects 
and V2 fluctuations, and such studies provide information on the magnitude of the systematic uncertainty. Non-flow 
effects are correlations not associated with the reaction plane, and include resonance decays. HBT correlations, flnal 
state interactions, and jets, to the extent that they do not participate in the flow. 

1. Event Plane method 

The essence of the Event Plane method [H is to first estimate the reaction plane. The estimated reaction plane is 
called the event plane and is determined by the anisotropic flow itself for each harmonic of the Fourier expansion of 
the anisotropic flow. The event flow vector Q2 and the event plane angle '^2 are defined by the following equations: 

Q2Cos(2*2) = Q2x = ^u;»cos(20,) (2) 

i 

Q2sin(2*2) = Q2y = ^WiSin(2(/),;) (3) 



= ( tan-1 ^ ) /2 (4) 



where the sum goes over all the particles i used in the event plane calculation. and Wi are the lab azimuthal angle 
and the weight for the particle i, respectively. In this analysis, the weights are taken to be the value of px in GeV/c 
up to 2 GeV/c and then constant at 2.0 above that px ■ 

The observed V2 is the second harmonic of the azimuthal distribution of particles with respect to this event plane: 



,obs 



(cos[2(0-vl/2)]) (5) 



where angle brackets denote an average over all particles with their azimuthal angle (/> in a given phase space. Since 
finite multiplicity limits the resolution in estimating the angle of the reaction plane, the real V2 has to be corrected 
for the event plane resolution by 

yObs 

" (cos[2(vI/2 - VI/,)]) 

where brackets denote an average over a large event sample, and vj/, is the angle of the reaction plane. The event 
plane resolution is estimated by the correlation of the event planes of two subevents. The event plane resolution for 
the subevents with the assumption of pure flow correlations between the subevents is 



(cos[2(vI/^ - *.)]) = y^(cos[2(f^-vI/f)]) (7) 

where A and B denote two subgroups of tracks. In this analysis, we use two random subevents with equal numbers 
of particles. Further, the full event plane resolution is obtained from the resolution of the subevents: 

(cos[2(*2 - = C(cos[2(*^ - *,)]) (8) 

where C is a constant calculated from the known multiplicity dependence of the resolution 0]. In the case of low 
resolution ( < 0.5), C is equal to y/2 0. The actual event plane resolutions, for the centrality bins used in this 
analysis of 0-10%, 10-40% and 40-80%, were 0.658 ± 0.0006, 0.818 ± 0.0002 and 0.694 ± 0.0004, respectively 
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2. ri-subevent method 



The 77-subevent method attempts to reduce the contribution from non-flow effects (mostly due to short range 
correlations) by correlating particles separated in pscudorapidity. This technique is similar to the event plane method, 
except one defines the event fiow vector for each particle based on particles measured in the opposite hemisphere in 
pseudorapidity: 

r . (COS[2(0^^ ~^2,„^)]) 

V2{v±} = I (9) 

^(cOs[2(vl/2,^^-Vl/2„J]) 

Here ^'2,77+ (^2,7)_) is the second harmonic event plane angle defined for particles with positive (negative) pseudora- 
pidity. An 77 gap of I 77 I < 0.075 between positive and negative pseudorapidity subevents is introduced in order to 
guarantee that non-flow effects are reduced by enlarging the separation between the correlated particles. In Eq. ^ 
the non-flow effects (correlations) are reduced in both the observed flow (numerator) and the event plane resolution 
(denominator). Depending on the nature of the remaining non-flow effects, V2 measured this way may have values 
which are cither lower or higher than those obtained with the standard method. 



3. ^-particle cumulant method 



A method to calculate V2 from true four-particle correlations was developed in Ref. [26[ and it has already been 
used by STAR It uses the cumulant relation 

C{4} EE (u„4U„,2<,3<,4) - 2 {unsK^^y = -«rl{4} , (10) 

where u„.j = e*"*^^ . The cumulant allows one to subtract the two-particle correlations, including two-particle non-flow, 
from the four-particle correlations. In practice, cumulants are calculated using the generating function from Ref. [26j 
and described in Ref. [2^ : 

where z = |z|e'" is an arbitrary complex number, with z* denoting its complex conjugate. The cumulants are related 
to the generating function by 

M . {{G^iz))'/'' - l) = E ||7^{2fc} . (12) 
The fit to the C{4} term is needed. The fourth root of the negative of it gives W2{4}. 



4- Lee-Yang Zero method 



The Lee- Yang Zero method [l^, [IHl is based on a 1952 proposal of Lee and Yang to detect a liquid-gas phase 
transition. As opposed to the four-particle cumulant method which is sensitive to the correlations of four particles, 
this method is sensitive to the correlations of all the particles. Thus it is supposed to remove non-flow correlations 
to all orders. It has so far been used only to analyze one set of experimental data (3^ and one set of transport 
calculations [33| . The method utilizes the second- harmonic flow vector, Q2, projected on to an arbitrary laboratory 
angle, 6: 

M 

Q^ = ^u;,cos(2(<^,-0)), (13) 
i=i 

where the sum is taken over all the particles j with lab angles 0j and weights Wj . For this method the weights are 
taken to be the value of pr in GeV/c for unidcntiflcd charged hadrons and 1.0 for identified particles. We have taken 
five equally spaced values of 9 to average out detector acceptance effects. The results were not different when 20 



8 



values of 9 were used. The theory of the method [24[ is to find a zero of a complex generating function, but in practice 
the first minimum of the modulus of the generating function along the imaginary axis is used. The sum generating 
function based on Q2 is given by 

G'Sr) = I (e^''^') I, (14) 

where r is a variable along the imaginary axis of the complex plane and the average is taken over all events. When 
data are analyzed in small batches, the G^^^iii") histograms arc combined before finding the first minimum. Such a 
histogram is shown in Fig. [D^a). The square of the modulus is used to determine the first minimum. The position 
along the imaginary axis of the first minimum of the modulus of the generating function at the lab angle 6 is called 
Tq, and is related to the "integrated" flow by 

V^^3,i/rl (15) 



V2 = {V^)JM, 



(16) 



where joi = 2.405 is the first root of the Bessel Function Jq and M is the multiplicity. In the second equation the 
average is taken over the lab angles 9. However, Eq. (fTB|) is only valid for unit weights. Normally, the anisotropic fiow 
parameter averaged over pT and rj is obtained by taking the yield-weighted average of the differential flow. For unit 
weights this has been shown to agree with the "integrated" flow from Eq. (fTB)) . The differential fiow obtained by a 
second pass through the data is given by 



cos[2m( 



^ Jm{j 01) 



Re 




(17) 



where m — \ for V2 and m = 2 for W4. The average in the numerator in the first equation is over the particles of 
interest and the average in the denominator is over all events. 





FIG. 2: (Color online) Examples of the modulus of the second harmonic Lee- Yang Zero Generating Functions plotted as a 
function of the imaginary axis coordinate, r. The Sum Generating Function is shown in (a) and the Product Generating 
Function in (b). The vertical arrows indicate the positions of the first minimum, called ro. Note that in (b) the horizontal scale 
does not go out as far because the calculations were terminated. All data are from ys„„= 200 GeV Au -\- Au collisions. 

The product generating function is 

M 

Gli^r) = I (Hfl + ^rw, cos(2(0, - 9))]) \ . (18) 

This takes more computer time because the product over all particles has to be calculated for each value of r. An 
example is shown in Fig.[2{b). It can be seen that the sum generating function oscillates after the first minimum, but 
the product generating function rises very fast. Thus, for the product generating function, the calculation was halted 
when |G'2(jr)p got larger than 1000. This happened at various r values between 0.2 and 0.4. While the method using 
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the sum generating function is slightly faster than the standard method 0, using the product generating function is 
about four times slower. For the product generating function, Eqs. p51 [TB]) still hold, but the differential flow is given 

by 

where again the average in the numerator is over the particles of interest and the average in the denominator is 
over all events. Although the sum generating function works fine for V2, analyses for (and vi) have to be based 
on the product generating function [25]. This is because the product generating function is better at suppressing 
autocorrelation effects which are more important for mixed harmonics. All methods used in this paper have been 
tested on simulated data. Also, since drift of the beam in the detector over time might simulate the effect of anisotropic 
flow, run-by-run recentering of the Q vector was applied, but produced no improvement in the results. 

The errors were calculated from the variation of the results for different event sub-samples. For very large errors 
this technique could underestimate the error because even when there is no flow the method will find a minimum from 
a fluctuation. In fact the Lee- Yang Zero method only works for sufficient signal-to-noise ratio. Since the signal is V2 
and the noise is proportional to 1/Vm, the parameter x ~ V2\/M determines the applicability of the method. We 
find that the errors get large and the results scatter when x < 0.8, and thus the results are presented here only for 
10-50% centrality. The method fails for more central collisions because V2 is small, and for more peripheral collisions 
because the multiplicity is small. 



III. RESULTS 



A. Charged hadrons 



To evaluate the different flow analysis methods and to estimate systematic uncertainties, charged hadrons were 
analyzed first. Figure [3] shows V2{r]) for both the Lee- Yang Zero and Event-Plane methods. For the Event-Plane 
method, the event plane was taken from the main TPC, both for tracks in the TPC as well as the FTPC. For the 
Lee- Yang Zero method, where there is no event plane, tracks in all three TPCs were used. The Lee- Yang Zero results 
arc for the product generating function but are in agreement with the sum generating function results. Elliptic fiow 
falls off in the Forward TPC (FTPC) covering | ry | from 2.6 to 4.2. This fall off is probably because the spectra 
as a function of pr are steeper at high -q and give less weight to the large V2 values at high pr [131 ■ Agreement 
in the FTPC region between two-particle and multiparticle methods has been seen previously [8|. Having a gap in 
pseudorapidity between the particles being correlated reduces the non-flow effects due to short range correlations. 
Indeed, PHOBOS correlates particles with an event plane from a different part of their detector, which is essentially 
similar to the ?7-subevent method. With | ?7 |< 1, PHOBOS [IHl data points are consistent with STAR Lee- Yang Zero 
data although it appears that PHOBOS data may be more peaked. Averaging over the TPC rj region | 77 | < 1.0, the 
V2{pt) values are shown in Fig. |4] (a) together with Event-Plane and 4-particle cumulant results. For these charged 
hadrons, the ratio of the four-particle cumulant result to the Event-Plane method shown in Fig. [3] (b) falls off as pt 
increases. This indicates a non-flow effect in the Event-Plane method which increases with pr as one would expect 
for the contribution of jets. On the other hand, the Lee- Yang Zero ratio seems to be flat. Figure [5] (a) shows the 
integrated V2 (averaged over rj and px) for four different analysis methods as a function of centrality. The ratios to 
the Event-Plane method are shown in Fig. O (b). The centrality region where all methods have reasonable error bars 
is 10 to 50%. The Event-Plane method appears to be about 15% higher compared to the other methods known to 
greatly reduce non-flow effects. This effect was already seen in the differential data as a function of rj for px < 2 
GeV/c in Fig. [3] and as a function of pr in Fig. [4l For the most peripheral collisions non-flow might be larger and 
for the most central collisions fluctuations could be important. From Fig. Owe would estimate the systematic errors 
at these other centralities to be 20%, and also probably this same value for minimum bias events. 
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FIG. 3: (Color online) V2 for charged hadrons from the Lee- Yang Zero Product Generating Function (solid circles) and from 
the Event-Plane method (open circles), as a function of pseudorapidity. Both sets of data have been averaged over pr from 
0.15 to 2.0 GeV/c and centrality from 10 to 40% of ^8^:^= 200 GeV Au + Au collisions. For comparison, the PHOBOS data 
(10-40%) [m (crosses) are also shown. The error bars are shown only for the statistical uncertainties. 
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FIG. 4: (Color online) (a) V2 as a function of pr for charged hadrons with | 77 | < 1.0 in 10-40% Au-'rAu collisions, at y/s^= 
200 GeV, from the Event-Plane method (open circles), 4-particle cumulant method (solid squares), and Lee- Yang Zero method 
(solid circles) with Sum Generating Function, (b) The ratios to the polynomial fit to U2{EP} are shown for i'2{4}/i'2{EP} and 
t'2{LYZ}/w2{EP} as a function of transverse momentum. The error bars are shown only for the statistical uncertainties. 
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FIG. 5: (Color online) (a) pT-integrated charged hadron V2 in the TPC as a function of geometrical cross section. Shown are the 
Event-Plane method (i;2{EP}) (open circles), Lee- Yang Zero method with Sum Generating Function (solid circles), Lee- Yang 
Zero method with Product Generating Function (open stars), and 4-particle cumulant method (u2{4}) (solid squares). For 
the TPC, I 77 I < 1.0 was used, except for Lee- Yang Zero method with Product Generating Function where the rj limit went to 
1.3. (b) V2 divided by i;2{EP}. All data are from y^s^— 200 GeV Au + Au collisions. The error bars are shown only for the 
statistical uncertainties. 
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B. K°s and A 

In order to estimate the particle dependence of systematic errors, mesons and A baryons were analyzed with 
different flow analysis methods. Figure [H] shows 10-40% V2{pt) of (a) A + A and (b) Kg obtained with the Event- 
Plane, Lee- Yang Zero, and 77-subevent methods. Ratios of V2 from these various methods to the Event-Plane method 
are shown in Fig. [S] (c) and (d) for A -I- A and , respectively. The results from Lee- Yang Zero method are 
about 10% lower than those from the Event-Plane method. This is consistent with charged particles, indicating 
that non-flow effects for and A are also reduced in the Lee- Yang Zero method. However, the ratio of the Lee- 
Yang Zero to Event-Plane method appears to be flat up to 5 GeV/c, considering the large statistical uncertainties. 
This is again similar to the trend observed for charged particle Lcc-Yang Zero results, which are shown as shaded 
bands. The results for the Ty-subevent method are higher than for the Event-Plane method especially at low px- 
The V2{r]}/v2{EP} ratio can be better understood by factorizing the ratio into an observed V2 term and a resolution 
term {V2'"' {rj} / V2^^ {F,P})x{resolution{F,P} /resolutionlrj}). In this analysis we calculate the 77 subevent resolution 
by correlating the event planes from the different 77 hemispheres. In this case, the non-flow effects are reduced in 
both V2'^ and the resolution. These factors contribute in the opposite direction to V2{r]}/v2{EP}; non-flow in the 
resolution term increases the ratio while non-flow in the ^2^^ term decreases the ratio. The V2{r]} / V2{EP} ratio is 
greater than unity because the resolution is more sensitive to non-flow than tig^'* and the decrease with px is caused 
by the increase of non-flow effects in W2''*{EP} with increasing px- 
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FIG. 6: (Color online) 112 as a function of pr for 10-40% centrality using Event-Plane method (open circles), Lee- Yang Zero 
method with Sum Generating Function (solid circles), and 77-subevent method (open crosses), are shown in (a) and (b) for A 
and Kg, respectively. All data are from y^s^= 200 GeV Au + Au collisions. The ratios, ?;2{LYZ}/u2{EP} and V2{ri}/v2{EP}, 
are shown in (c) and (d) for A and Kg, respectively. D2{LYZ}/d2{EP} ratios for charged hadrons are shown as shaded bands. 
The error bars are shown only for the statistical uncertainties. 



C. Systematic uncertainty of A from feed-down 

We estimate the reaction plane orientation from the azimuthal distribution of charged particles measured with the 
TPC (I77I < 1.3), constructing the second harmonic event plane flow vector Q2. The TPC is also used to reconstruct 
A and A liyperons via their charged decay daughters, tt^ or p{p). Elliptic flow of A and A is measured by correlating 
the hyperon azimuthal angle with the event plane vector Q2. Although the correlation strength is mainly defined by 

the hyperon elliptic flow, w^'^, in such an approach other sources of correlations (non-flow effects) may contribute 
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and bias the measured V2'^ values. 

We do not distinguish between A and A particles produced from the secondary decays (for example, S~ — > A + 7r~ 

(S^ ^ A + 7r+), S](1385)" A + tt" (S(1385)+ ^ A + 7r+) or 1]° A + 7 ) and hyperons which originate directly 
from the primary interaction. Indirect hyperons lead to the presence of extra correlations that are not related to 
the reaction plane between hyperons and other charged particles produced in the collision. Note that the charge 
combinations for these correlations are opposite for A and A particles. 

To estimate the contribution of these non-flow correlations from hyperon feed-down effects we use the charge 
subevent technique. For this method we introduce two event plane vectors: Q2 constructed from positively charged 
particles and Q2 from negatively charged particles. We then estimate the contribution to the A and A elliptic flow 
measured with the full event plane vector Q2 by considering the following ratio: 



^i?{FeedDown} = ^im±iMl. (20) 
v^{Qj}+v^{Q2-} 

Here V2'^{Qf} denotes A (A) elliptic flow values measured from correlations with Q^. The numerator in Eq. pO]) 
contains the contributions from non-flow correlations attributed to feed-down effects, while the denominator is free of 
this. From this study we found that contribution of non-flow effects from feed-down of secondary A and A hyperons 
is < 2%. 



D. Other systematic uncertainties 



In order to estimate the systematic uncertainties in the identified hadron V2, we employed the standard Event- 
Plane method, the 77-subcvent method, and the Lee- Yang Zero method. The results of the analysis for A-l-A and 
Kg are shown in Fig. [B] The ratios to the event plane results are shown in the lower panels. We limited ourselves 
to the 10-40% centrality bin where the results from the Lee- Yang Zero method are most reliable. For comparison, 
the charged hadron results from the W2{LYZ}/'i;2{EP} ratios (Fig. |4]) are also shown in the figure as shaded bands. 
There is no clear pT dependence of the ratio ■i;2{LYZ}/u2{EP} for either A+A in (c) or for Kg in (d) although, within 
the statistical errors, similar trends are seen for both identified particles. The overall systematic errors are on the 
order of 15% in the pt region studied. The subevent method, however, introduced an additional factor described in 
Section Hi C 21 that leads to the enhanced ratio for both A+A and Kg. This opposite effect is also within the order of 
15%. 

In order to obtain good statistics, the Event-Plane method is used for most of analyses of identified particles. 
Depending on the analysis method, systematic uncertainties from variations in particle identification cuts, background 
subtractions, and summing of centrality bins are also estimated. By varying particle identification (PID) cuts, which 
change signal over background ratios by a factor of 3, the systematic uncertainty from the PID cuts is estimated to 
be about 5% below 4 GeV/c. From 4 GeV/c to 6 GeV/c, this effect is larger in central collisions than peripheral 
collisions, and in the 0-10% bin it is about a 10% effect. To estimate the uncertainty from background subtractions, 
background variations from different second and fourth order polynomial fit functions in Fig. [T] were propagated to 
measured V2 values with Eq.[TJ The effect is less than 3%. When combining ccntralities, the combined V2 value should 
be a yield- weighted average of V2 values in small centrality bins. In the method used previously 0, f2 values are 
taken as weighted observed V2 corrected by weighted event plane resolution. With the V2 versus rriinv method, similar 
corrections were calculated. They are less than 5% below 6 GeV/c. 

In summary, for charged particles, a 15% difference for {V2) at mid-rapidity for 10-40% collisions between the 
event plane and the Lee- Yang Zero methods has been observed. The difference between ti2{EP} and V2{4} is smaller, 
~10%, but for the more peripheral and more central collisions it seems to be closer to 20%. For V2{pt) of Kc^ and A, a 
difference also is observed between the Event-Plane and the Lee- Yang Zero methods. However, the comparison with 
charged particles shows that, within the much larger statistical uncertainties for the A and Kg analysis, the different 
magnitudes of the estimated flow from the different methods observed are similar for charged particles and for the A 
and Kg particles. The uncertainty used is 15%. In a latter section, the results are presented for S~ and the fJ^ . For 
these particles, due to limited statistics, only the Event-Plane method has been used and therefore no real estimate of 
the systematic error is available. Instead what was done, was to show the estimated systematic uncertainties obtained 
for the charged particles, taking into account effects from background estimation, summing centralities, and variations 
in cuts. A summary of our best knowledge of systematic errors is given in Table [D 
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centrality 


0-80% 


40-80% 


10-40% 


0-10% 


charged hadrons 


20% 


20% 


15% 


20% 


identified particles 


N/A 


N/A 


15% 


N/A 



TABLE I: Systematic errors of V2 estimated from the difi'erent flow methods summarized as a function of centraUty for y^s^ = 
200 GeV Au + Au collisions for unidentified charged hadrons and identified particles. The px region covered for charged hadrons 
is 0.5 < pt < 7.0 GeV/c and for identified particles Kg, p, and A is 0.5 < pr < 5.0 GeV/c. 



IV. DISCUSSION 
A. Charged hadrons 

Figure [5] shows that for the Lee- Yang Zero method the sum and product generating functions agree, but are shghtly 
lower than the four-particle cummulant method. The Event-Plane method appears to be about 15% higher than the 
other methods. This could be due to either non-flow contributions increasing the event plane results or fluctuations of 
V2 decreasing the multi-particle methods Q . Charged hadron results give an indication of the systematic uncertainty 
inherent in the Event-Plane method, which is used for most of the identified particles in this paper in order to reduce 
statistical errors. Thepr dependence of this effect can be seen in Fig. 01 Note that the mentioned 15% non-flow effect 
is extracted only from the 10-40% centrality window. In other centrality bins, the effect may be larger. It should be 
possible to study this effect as a function of pT for all centralities with a future higher statistics data sample. 



B. Identified hadrons 

The results for 7r+4-7r~, p + p, Kg, A+A , and il^+il^ are shown in Figure [7] for various centralities of 

Au + Au collisions at 200 GeV. Shown are results for minimum bias and three other centrality bins. All V2{pt) results 
are from the Event-Plane method. The systematic uncertainties extracted from PID cuts, background subtractions, 
and combining centralities, are shown as shaded bars in the figure. The systematic uncertainty in the method itself 
is not included. The shaded band in plot (c) indicates the systematic uncertainties for and A for the 10-40% 
centrality bin, as discussed in the previous section. In addition, the systematic uncertainties for pions and protons 
from Event-Plane and 77- sub event analysis are plotted as the shaded band in plot (a). The results from an ideal 
hydrodynamic model [Sy, [s^l are displayed by the lines. 

Figure [7] shows that the ideal hydrodynamic model calculations reproduce the mass ordering of V2 in the relatively 
low Pt region (the heavier the mass, the smaller the V2) but overshoot the values of V2 for all centrality bins. There 
seems to be a pT dependence in the disagreement, and for more central collisions, the overshoot does not take 
place until a higher px- In other words, the system agrees better with the ideal hydrodynamic model for more 
central collisions. Although we do not expect a large non-flow contribution at the low transverse momentum region, 
the centrality selections between the model calculations based on the impact parameter and the data based on the 
multiplicity are different, which may also affect the model and data agreement. Note that we observe possible negative 
values of V2{pt) for the heavier hadrons at the lowest observed pt in the most central Au + Au collisions. 

At higher px, the hydrodynamic type mass ordering evolves into a hadron type ordering (baryons versus mesons). 
There the results show two groups depending on the number of quarks in the hadron; the baryons are higher than the 
mesons. The effects are clearly shown in plots (a), (b) and (c). In the most central bin (d), however, the effect is less 
pronounced. For all pt, V2 evolves toward larger values in going from central collisions to more peripheral collisions. 
The ideal hydrodynamic model also predicts this centrality dependence though it fails to describe the behavior at 
higher pT- 

Figure [S] shows the same results as in Fig. [7] but as a function of the transverse kinetic energy Kt = rnx — m = 
■\/p|7+m^ — m. Here m is the particle mass. In this case, V2 for all hadrons at low Kt follow a universal curve, 
which appears to be monotonically increasing and almost linear in all centrality bins. The observed increase is slowest 
for the most central 0-10% bin. The corresponding results from the ideal hydrodynamic model calculations are also 
shown in the figure. The mass ordering in the model calculation is reversed when one plots V2 versus mx — m: the 
higher the mass the larger the value of V2. While the data seem to show a scaling in the low tot — w region, the 
model results do not show any scaling. 

In Figs.[9l[T0l and [TT] we discuss the properties of the centrality dependence of the observed scaling including both 
the low Pt < 2 GeV/c and the intermediate 2 < < 5 GeV/c regions. Figure [9] shows V2 scaled by the number of 
constituent quarks, V2lnq, for all strange hadrons including the pure multi-strange hadrons and 17. The left panels 
show the results as a function of pt scaled by the number of quarks, pT/riq, and the right panels as a function of 
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FIG. 7: (Color online) V2 of (open circles), A (open squares), H (filled triangles) and Q, (filled circles) as a function of pr for 
(a) 0-80%, (b) 40-80%, (c) 10-40% and (d) 0-10% in Au + Au collisions at ^/s^= 200 GeV. The error bars represent statistical 
uncertainties. The bands on the data points represent systematic uncertainties as discussed in the text. For comparison, pion 
(stars) and proton (filled squares) results are shown in (a). The systematic uncertainty of non- flow for K°c; and A for 10-40% 
(c) is plotted as a shaded band near 0. For comparison, results from ideal hydrodynamic calculations [36l. |37| are shown: at a 
given pt, from top to bottom, the lines represent the results for n, K, p, A, H, and fi. 

Kx/nq. Plots (a) and (b) are the corresponding scaled results. It appears that the scaling works better when the 
data are plotted as a function of transverse kinetic energy Kt, as in plot (b). The ideal hydrodynamic results are 
also shown in both presentations. Clearly, the hydrodynamic distributions are also better scaled when plotted versus 
Kt- Polynomial fits were made for all hadrons. The results arc shown as dot-dot-dashed lines in plot (a) and (b). 
The ratios of the data and the hydrodynamic lines over the polynomial fit arc shown in plots (c) and (d). For all 
data, there is scaling at pr/nq > 0.7 GcV/c or Kr/nq > 0.2 GcV/c^. The errors from the multi-strange hadrons 
(i hj and arc large, sec plots (c) and (f), but arc consistent with the scaling. This observed riq-scaling provides 
strong evidence that these hadrons are formed via a coalescence process at the end of the partonic evolution [ll], [I!] . 
Comparing to the light non-strange hadrons, strange hadrons participate much less in later stage hadronic rescattering 
processes |3g; thus, these distributions directly reflect the early dynamics of the collision at RHIC. It is interesting to 
note that the ideal hydrodynamic results scale neither at low px nor intermediate px- Therefore the observed scaling 
cannot be a general characteristics of hydrodynamic model calculations [l7j , although such calculations do show the 
observed mass ordering in the low pt region. 

Figures [TOl and [Til show the centrality dependence of the scaling properties and the ratios, respectively. Similar to 
the observations from Fig. [SI the conclusions from Figs. [TUI and [TT] are: 

1. There is a clear numbcr-of-quark scaling at intermediate px and better scaling in Kt for all hadrons studied 
here, but no scaling is observed at low pT- 

2. The ideal hydrodynamic model results do not show any scaling over the region 0.2 < pt < 5 GeV/c. 

3. These results are true for all centrality bins. 
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FIG. 8: (Color online) V2 in Fig. [7] is re-plotted as a function of rriT — m for K% (open circles), A (open squares), S (filled 
triangles) and (filled circles), for (a) 0-80%, (b) 40-80%, (c) 10-40% and (d) 0-10% in Au + Au collisions at ,Js;^= 200 
GeV. For comparison, pion (stars) and proton (filled squares) results are shown in (a). The ideal hydrodynamic calculations 
in Fig. [7] are also shown as a function of rtiT — m. 
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FIG. 9: (Color online) Number-of-quark scaled V2 {v2/nq) of identified particles versus pr/nq in (GeV/c) (left column) and 
{rriT — m)/nq in (GeV/c^) (right column). Dot-dot-dashed lines are the results of 6th order polynomial fits to J^g, A, H, and 
n. The ratios of the data points over the fit are shown in panels (c) and (d) for K^, tt, p, A, and H, and in panels (e) and (f) 
for n and cj) [3- The error bars are shown only for the statistical uncertainties. Ideal hydrodynamic calculations for vr, 
A, H, and n are presented by solid lines, dashed lines, dot-dashed lines, dot-long-dashed lines, and dot-dot-dot-dashed lines, 
respectively. The ratio of hydrodynamic calculations over the fit are also shown for comparison in panels (c) and (d). The data 
are from minimum bias (0-80%) Au -\- Au collisions at y^s^ = 200 GeV. 
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FIG. 10: (Color online) Centrality dependence of the number-of-quark scaled V2 {v2/nq) of identified particles versus pr/nq 
in (GeV/c) (left column) and {niT — m)/nq in (GeV/c^) (right column). The error bars are shown only for the statistical 
uncertainties. The 6th order polynomial fits are shown as dot-dot-dashed-lines. Ideal hydrodynamic curves [s^l are also 
plotted. All data are from ^/s^= 200 GeV Au + Au collisions. 
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FIG. 11: (Color online) Centrality dependence of the V2/nq ratio to a common polynomial fit is shown versus pr/riq in (GeV/c) 
(left column) and (mr — m)/nq in (GeV/c^) (right column). The error bars are shown only for the statistical uncertainties. 
Ideal hydrodynamic calculations over the same fit are also shown for comparison. Ideal hydrodynamic calculations for tt, K, 
A, H and Q are presented by solid lines, dashed lines, dot-dashed lines, dot-long-dashed lines, and dot-dot-dot-dashed lines, 
respectively. All data are from a/s^= 200 GeV Au + Au collisions. 
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(^trig 1 0'geom 


0-80% 


40-80% 


10-40% 


0-10% 


62.4 GeV 


£part 
^ part 
^hin 


0.3919±0.0003 
122±3 
249±13 


0.5426±0.0004 

39±5 
50±10 


0.2927±0.0003 
167±7 
338±18 


0.1108 ±0.0002 

320±3 
797±9 


200 GeV 


£part 
^ part 
Nbin 


0.3843±0.0001 
126±8 
293±36 


0.5343±0.0002 

42±7 
57 ± 14 


0.2829±0.0001 
173±10 
393± 47 


0.1054±0.0001 
326 ± 6 
939 ± 72 



TABLE II: List of participant eccentricity Spart, number of participants Npart, and number of binary collisions Nbin, from a 
Glauber calculation [39l . |40| ] for minimum bias and three other centrality bins. The errors are statistical from the calculations 
only. The biggest systematic error is probably from the centrality binning based on the impact parameter compared to the 
binning of the data based on the multiplicity. All parameters for 62.4 GeV Au+Au collisions are calculated in a similar fashion 
as for 200 GeV collisions. The p+p cross sections used were 36 mb at 62.4 GeV and 42 mb at 200 GeV. 



C. Universal scaling? 

In order to analyze the centrality dependence of the scaling properties, we normalize the riq-scaled elliptic flow V2 
by the participant eccentricity Spart from a Monte Carlo Glauber calculation [s^, |43|. (See Table [III for £part ) The 
results are depicted in Fig. [T^l The plots (a) and (b) show the doubly scaled quantities from three centrality bins as 
a function of px/nq and (tot ~ fTi)l'>^q, respectively. Both plots show an initial rise and a turn over to a flat region in 
the higher pT region. It is interesting to see in (a) and (b) that at a given ceirtrality, the elliptic flow of all hadrons are 
scaled as observed in the minimum bias case (Fig. . After the geometric effect has been removed by dividing by Spart 
in (a) and (b), the build up of stronger collective motion in more central collisions becomes obvious in the measured 
elliptic flow. This is consistent with the ideal hydrodynamic model calculations, shown as lines in (b), although the 
model results are much closer together. However, clearly there is no scaling amongst different collision ccntralities. 
Neither our data nor the model results indicate universal scaling with eccentricity. A careful inspection of the results 
presented in Ref. [4l[ shows there is no disagreement between data; rather, the statement of the universal scaling in 
Ref. [4l[ is not supported by the data. 

To further clarify the issue, instead of dividing the measured V2 by the corresponding eccentricity Spart, we plot 
V2{mT — itl)/ {nq x {V2)) for Kg , A , and S in Fig. [T^] (c). The values of (W2) (see Table UlI)) are obtained by averaging 
V2 as a function of transverse momentum weighted with the measured spectra. As one can see in the figure, for a given 
hadron, this scaling seems to work better. However, different hadrons seem to have different values of V2, especially 
for the top 10% centrality bin at the higher to^- 

Figure [m (d) shows the doubly scaled V2 again. But this time, the integrated values of V2 are extracted from the 
measurements of unidentified charged hadrons {v2)ch at the corresponding centrality bins. In this case, it appears 
that the scaling works better. It is interesting to point out that at the most central bin, see inset in Fig.[T2](d), the 
values of V2 become negative at low px for all hadrons. This is most likely caused by the strong radial flow developed 
in central Au+Au collisions [i^. Similar behavior has also been observed in V2 of A at SPS [isj . 

D. Integrated V2/tpart versus collision centrality 



^trig / geom 


40-80% 


10-40% 


0-10% 




0.0735 ± 0.000163 


0.0576 ± 0.000064 


0.0283 ± 0.000112 




0.0707 ± 0.0008 ± 0.0013 


0.0513 ± 0.0005 ± 0.0013 


0.0212 ± 0.0011 ± 0.0012 


<!> 


0.0851 ± 0.0111 ± 0.0020 


0.0658 ± 0.0082 ± 0.0016 


0.0210 ± 0.0116 ± 0.0050 (0-5%) 


A 


0.0899 ± 0.0010 ± 0.0013 


0.0609 ± 0.0006 ± 0.0019 


0.0221 ± 0.0012 ± 0.0029 




0.0858 ± 0.0045 ± 0.0000 


0.0577 ± 0.0032 ± 0.0023 


0.0220 ± 0.0028 ± 0.0015 



TABLE III: The pT-averaged 112 of identified particles (particle + anti-particle) from three centrality bins in ,ys^^— 200 GeV 
Au+Au collisions. The Event-Plane method was used to extract the values of V2. Statistical and systematic errors are shown 
as the first and second errors, respectively. 
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FIG. 12: (Color online) 112 scaled by the number of quarks (n,) and participant eccentricity {spart), {V2/ [uq x Spart)), of 
identified particles (particle + anti-particle) versus (a) the scaled pr/uq and (b) (ttit — m)/nq for three centrality bins. For 
comparison, ideal hydrodynamic model calculations [37| are shown as lines in (b). In (c) is shown the data from (b) scaled by 
the integrated V2 of each particle, instead of Spart- In (d) is shown the data from (b) scaled by the integrated V2 of all charged 
hadrons. The insert in (d) expands the low rriT region. The error bars are shown only for the statistical uncertainties. All data 



are from ,/s7 



200 GeV Au + Au collisions. 



The integrated elliptic flow values in Table Hill were obtained from the measured V2{pt) and separately parameterized 
Pt spectra. The V2{pt) were integrated over px weighted with the yield distribution from functions fitted to the 
spectra. To extend V2 to low pT, a sixth order polynomial and a rig-inspired function [31 were used to fit V2- The V2 
values in the table are the average values from these two sets of parameterizations. The systematic uncertainties are 
taken as half of the differences between values from two sets of fits. The statistical errors as a function of pt are fitted 
with third order polynomials and folded with the yield distributions into the errors of the integrated V2- The spectra 
for and A are from Ref . [l^ and the spectra for S are from Ref. [iJl ■ Data for ^-mesons are from Ref. [l^ . 

The centrality dependence of the ratio of the integrated elliptic flow (Table HIT)) over the eccentricity {v2/Spart) 
for charged hadrons, , 0-meson [l^, A and S are shown in Fig. [T31 All these results are from the Event-Plane 
method and the number of participants is the average in the centrality bin. For comparison, results from an ideal 
hydrodynamic calculation [33,[iB| are also shown as dashed lines. This ratio, to some extent reflects the strength of the 
collective expansion. At more central collision, one would expect a stronger expansion, hence the larger value of the 
ratio. This is what one sees in Fig. [T3] for charged hadrons. For identified hadrons, the increasing trend as a function 
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FIG. 13: (Color online) Centrality dependence of v^jepart versus number of participants for charged hadrons (crosses), K% 
(circles), (\> (stars) [31, A + A (squares) and H + H (triangles). Both unidentified charged hadron and identified hadron 
were analyzed with the standard Event-Plane method. All data are from y^s^= 200 GeV Au + Au collisions. The data points 
are displaced slightly horizontally for clarity. The statistical uncertainties and the systematic uncertainties are shown as bars 
and brackets, respectively. Ideal hydrodynamic model calculations are also shown as dashed lines [stI . |46| for, from top to 
bottom, n, H, A, p, K, and tt. 



of Npart is there despite the large error bars. In the ideal hydrodynamic calculations [37|, |46[, the first order phase 
transition and freeze-out temperatures are set to be 165 MeV and 130 MeV, respectively. With these parameters, 
the ideal hydrodynamical model results describe the pion, kaon and proton transverse momentum spectra [stI . [46| . 
In a pure hydrodynamic model, one deals with energy-momentum cells rather than any specific type of hadrons, thus 
the initial condition, the equation of state and the freeze-out conditions used in the calculation arc the same for all 
hadrons. Such assumptions may not be applicable to all hadrons since some of them vifill continue to interact even 
after hydrodynamic freeze-out [431. 

As expected in an equilibrium scenario, the model results show little sensitivity to the collision centrality. However, 
it is interesting to note that there is a clear hadron mass dependence of vi normalized by £part from the model 
calculations which is not seen in the data. It is not clear whether the mass dependence is from the collective motion 
at early time or is the effect of the hadronization process in the calculation. On the data side, the errors are too 
large to allow comparisons with model results. As one can see in Fig. [131 after N^art ^ 170, the measured ratios for 
the strange particles approach that from the ideal hydrodynamic model calculations. The consistency between model 
results and data indicates that the system created in 200 GeV Au + Au collisions may reach local thermalization in 
central collisions when the number of participants is larger than ~ 170. 



E. Energy dependence 

The transverse momentum dependence of eccentricity-scaled vi ratios are shown in Fi g. 1141 Plots (a) and (b) show 
the ratio of data from 62.4 GeV [Au -I- Au collisions) and 17.2 GeV (Ph -I- Ph collisions |^|T over 200 GeV Au -f Au 
collision data, respectively. It has been observed [48| that since the charged multiplicity production per participant is 
proportional to the square root of the CM energy [4g, dN / d-q/ {Npart/'^) oc y^SNN, a stronger flow is expected from 
higher energy collisions [i^. In Fig. [TH the ratios for both Kg and A are similar. While the higher energy ratios 
show a decreasing trend as a function px, the lower energy ratios seem to increase with transverse momentum. In 
the low pt region (< 1.5 GeV/c), the strength of flow is similar in 200 GeV and 62.4 GeV Au + Au collisions. For 
the 0-80% Au-|-Au collisions, the values of the participant eccentricity from 62.4 and 200 GeV arc 0.392 and 0.384, 
respectively. The lack of energy dependence in V2 for Kg and A is due to the similarity in t he p articipant eccentricity. 
The PHOBOS experiment reported a similar observation for the V2 of charged hadrons [501. discussed in [48j . 
from 62 GeV to 200 GeV, the participant-normalized charged hadron density at mid-rapidity increased by about 50%. 
However, we do not observe a change of a similar size in the participant eccentricity and indicating that a large 
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(GeV/c) 

FIG. 14: (Color online) The pT dependence of the eccentricity-scaled V2 ratios of (a) 62.4 GeV and (b) 17.2 GeV [l^] over 200 
GeV data for Kg (circles) and A + A (squares). As indicated in (b), the SPS data points are from 5-25% and the RHIC data 
points are from 10-40%. Error bars are statistical only. 



fraction of the particle production occurs at the later stage of heavy ion collisions at these beam energies. 

The ratio is less than unity in the low px region (see Fig. [M] (b)) indicating that the flow is weaker in the lower 
energy Pb + Pb collisions [43]. The collective velocity parameters, extracted from the transverse momentum spectra, 
are also found to be larger in 200 GeV Au + Au collisions than those from 17.2 GeV Pb + Pb collisions [ll|. Since 
elliptic flow develops at a relatively early stage of the collisions, the observed increase in collective flow in Au + Au 
collisions at RHIC is therefore caused by early partonic interactions 0, H, Q ■ 

Figure [TSl shows the centrality dependence oiv2 normalized by the number of quarks and eccentricity {V2/ {nq£part)) 
for identified hadrons from (a) 62.4 GeV [5l| and (b) 200 GeV Au + Au collisions. Within error bars, data from both 
energies are similar. At the low transverse energy region, the scaled V2 shows almost a linear increase and then becomes 
fiat. For more central collisions the turning point is at higher values of [niT — m) / Uq. Recently, PHENIX has reported 
a charged hadron scaling with eccentricity, system size, and the transverse energy {jtit — m) up to 1 GcV/c^ [4l] |. 
As one can see from the figure, at a given centrality, independent of the collision energy, there is a clear scaling: all 
values of V2/{nqSpart) coalesce into a single distribution. On the other hand, it is clear in the figure that at different 
centralities the shape of the distributions are different, meaning that there is no scaling in the measured V2 with the 
eccentricity, especially in the higher transverse energy region. 



V. SUMMARY 



We present STAR results on the elliptic flow V2 of unidentified charged hadrons, strange and multi-strange hadrons 
from y/s^ = 200 GeV Au+Au collisions at RHIC. The centrality dependence of V2 over a broad transverse momentum 
range is presented. Comparison of different analysis methods are made in order to estimate systematic uncertainties. 
The rapidity dependence of the charged hadron V2 from these measurements is consistent, at both mid-pseudorapidity 
and forward rapidities, with both STAR and PHOBOS [s^ reported results. In particular, the results for V2 from 
the Lee- Yang Zero method for charged hadrons, , and A are shown for the flrst time at RHIC. The non-flow 
effects, studied in the 10-40% centrality window, are on the order of 10% within 0.2 < pr < 3 GcV/c and up to 25% 
at pt ~ 6 GeV/c. 

In the relatively low px region, pt l£ ^2, GeV/c, a scaling with niT — m is observed for idcntifled hadrons under 
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FIG. 15: (Color online) Centrality dependence of V2 / {nqe-part) for identified hadrons (particle + anti-particle) from (a) 62.4 
GeV [5l[ and (b) 200 GeV Au + Au collisions. Error bars are statistical only. For the 10-40% centrality the systematic 
uncertainty should be the same as in Fig. [T] 



study in each centrality bin and there is a clear centrality dependence in the scaling. However, we do not observe 
V2{mT — m) scaled by the participant eccentricity to be independent of centrality. The largest values of the participant 
eccentricity scaled V2 are in the most central collisions. For the most central collisions (0-10%), negative values of 
V2 at the lowest pt studied have been observed for both Kg and A . This is the first time a negative V2 in Au + Au 
collisions at RHIC has been found. It is consistent with the strong expansion observed in hadron spectra analysis [ll| . 
In the higher px region, 2 < px < 6 GeV/c, number-of-quark scaling is observed for all particles under study. For 
the multi-strange hadron f2, which does not suffer appreciable hadronic interactions, the values of V2 are consistent 
with both mx — m scaling at low px and number-of-quark scaling at intermediate px- 

As a function of collision centrality, an increase of pT-integrated V2 scaled by the participant eccentricity has been 
observed, indicating stronger collective flow in more central Au + Au collisions. However, in the higher transverse 
energy region there is no scaling of V2 with eccentricity. 

The energy dependence of V2 for and A was presented for y^s^ = 17 GeV Pb + Pb, and 62 and 200 GeV 
Au + Au collisions. No clear systematic trend was observed, but the differences were only on the order of 10%. 

For comparison, results from ideal hydrodynamical model calculations were used. The calculations over-predict the 
data at > 2 GeV/c. At low px , although the model predicts correctly the mass hierarchy observed in the data, 
there is no scaling with mx — m, and no scaling with the number of quarks is observed throughout the px region for 
all hadrons. We observe that the mass ordering at low px alone is not sufficient to claim thermalization in Au + Au 
collisions at RHIC. 
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